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SUMMARY 

A  method  of  obtaining  the  bivariate  cumulants  of  any  order  is 
given  for  a  truncated  bivariate  normal  distribution  where  one  of  the 
variates  is  truncated  at  w^.  Some  representative  values  are  displayed 
in  tables .  , ,  C  >  W  j"  . 

Some  key  words:  Bivariate  cumulant;  Singly  truncated  bivariate 
normal  distribution. 


1 .  INTRODUCTION 

Let  the  joint  density  of  a  standardized  bivariate  normal  distri¬ 
bution  be  given  by 

$(x,y;p)  *  (2ir)_1(l-p2)  ^  exp{ -± (x2  -  2pxy  +  y2)/(l-p2)} 
for  -  «  <  x,y  <+  co  , 

We  will  also  use  the  notation 

G’(x)  •  (2ir)  ^  exp{-£x2}, 
and  G(x)  *  fX  G’(t)dt 

'  —  09 

for  -»  <  x  <  +  «  , 

for  the  standardized  univariate  density  and  cumulative,  respectively  . 
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Then  if  the  X  variate  is  truncated  below  wq,  the  joint  density  of  the 
singly  truncated  bivariate  normal  distribution  (STBVND)  is  given  by 

f(x,y;p)  *  ij>(x,y;p)/G(-wo)  for  wq<x<+“,-“<  y  <  +  «®. 
Our  purpose  is  to  obtain  the  bivariate  cumulants  corresponding  to 
f (x,y;p) . 


2 .  BIVARIATE  CUMULANTS 

Cook  (1951)  illustrated  three  methods  of  deriving  bivariate  cumu¬ 
lants.  Cumulants  of  all  orders  of  the  bivariate  distribution  may  be 
worked  out  by  choosing  the  appropriate  operation.  She  gave  all  the 
formulae  for  bivariate  cumulants,  up  to  i  +  j  =  6,  As  the  order 

of  the  cumulants  increases, the  number  of  terms  increases  greatly. 
Johnson  and  Kotz  (1972)  gave  bivariate  cumulants  only  up  to  i  +  j  *  2. 
Gajjar  and  Subrahmaniam  (1978)  obtained  bivariate  moments  up  to  order 
4.  We  give  here  a  general  formula  for  bivariate  cumulants  for  any 
order . 


The  moment  generating  function  of  a  STBVND  is  given  by 
M(tl’C2^  *  [G(-wo)]-1G(t1  +  pt2  -wq)  exp  [r(t^  +  2pt1t2  +  t2)  ] . 
where  t^  corresponds  to  x  and  t2  corresponds  to  y.  Since  the  cumulant 
generating  function  is  ln[M(t^ , t2) ] ,  we  have  the  following  expression 
for  the  cumulant  generating  function  K(t^,t2) 

K(t1,t2)-- ln[G(-wQ)]  +  ln[G(t1+pt  -w  )]  +  t*/2  +  pt^  +  t^/2 

and  the  cumulant  is  obtained  by  taking  the  i-th  partial  derivative 

with  respect  to  t^  and  the  j-th  partial  derivative  with  respect  to 


t2  and  setting  t^ 


0. 
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We  obtain 

k  -  G'(w  )/[G(-w  )] 
1U  o  o 


*01  "  pK10 


Sn  *  w  G*(w  )/[G(-w  )]  -  [G’(w  )/G C-w  )T+1 
ZU  O  O  o  O  o 


=  VlO  ~  *10  +  1 


k02  =  p  ^(f15  +  1 


and 


*ij  =  P 


3  3 


1+j  -1 


(g* (x)1 


i+j-1 


3x 


G(x) 


for  i  +  j  i  2 


X  =*  -w 


The  expressions  for 


1  fG’ (x)  1 


3x 


G(x) 


become  very  cumbersome  as  is 


illustrated  by  the  following. 
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[g'  (x)1 

fG'  (x)^ 

[g’(x)12 

3x 

[g(x)  j 

a  ”*X 

[g(x)  j 

[G(x)  J  ’ 

3 

ax2 

ai 

3x3 


G'(x) 


(G(x)  J 
G 1  (x) 


(x2-l) 


fG’(x)1 

+  3x 

fG'  (x)l 

2 

O 

fG’(x)] 

3 

lG(x)  J 

,G(x)  j 

T  Z 

[G(x)  j 

* 

G(x) 


(-x  +  3x) 


-  12x 


rG'(x) 


|G(x)  J 


-  (7x  -  4) 

4 


fG r (x)l 2 


G(x) 


[G’(x)l 

>  3 

—  & 

fc’(x)] 

(g(x)  j 

J  “  0 

,G(x)  , 

Hence,  we  look  for  a  way  to  generate  these  derivatives  recursively. 
We  let  g  *  g(x)  -  G'(x)/[G(x)]  and  h  -  1/g  -  G(x)/[G'(x)]  and  note 
that  the  derivatives  of  h  are  easily  obtained  as 
h'  -  1  +  xh 

h(n)  -  (n-l)h(n~2)  +  xh(n"1)  for  n  >  2,  where  h(0)  =  h. 


4 


Since  gh  *  1,  an  application  of  Liebniz's  rule  gives 

I  A  -  o 

j-0  J 


or 


8(n>  -  -8  l  A  8(”'3)hCJ). 

i-1  3 


Hence,  we  can  obtain  any  derivative  of  G'(x)/G(x)  in  terms  of. lower 
order  derivatives  and  derivatives  of  h  which  are  given  by  the  recur¬ 
sion  formula  for  h^  . 


3 .  APPLICATION 

In  the  problem  of  screening  based  on  a  singly  truncated  bivariate 
normal  distribution,  one  needs  to  know  the  distribution  of  the  sample 
correlation  coefficient.  This  can  be  achieved  by  supplying  the 
bivariate  cumulants  given  above  to  Gayen's  (1951)  results. 


4 .  TABLE 

In  the  accompanying  table  we  give  values  of  k^q,  *30’  <40 

and  These  cumulants  are  independent  of  p.  However,  the  remaining 

cumulants  up  to  order  5  may  be  obtained  from  these,  using  the  following 
simple  formulas  which  were  obtained  from  the  general  formula  for 
given  above: 

*01  "  p*10 

*02  *  p  ^20  ^  +  1 


*11  “  p  *20 

Kij  •  pJ'l+j,0  for  1  +  3  ”  2 


TABLE  1 


BIVARIATE  CUMULANTS  OF  A  SINGLY  TRUNCATED 
BIVARIATE  NORMAL  DISTRIBUTION 


W 

o 

*10 

*20 

*30 

K,  . 

40 

*50 

-3.0 

.004438 

.986667 

.035680 

-.081046 

.139672 

-2.8 

.007936 

.977717 

.054810 

-.110766 

.154800 

-2.6 

.013647 

.964333 

.080062 

-.141556 

.148857 

-2.4 

.022580 

.945299 

.111173 

-.168418 

.114778 

-2.2 

.035975 

.919561 

.146778 

-.185535 

.052033 

-2.0 

.055248 

.886452 

.184395 

-.187855 

-.031092 

-1.8 

.081893 

.845887 

.220752 

-.172786 

-.118818 

-1.6 

.117352 

.798466 

.252404 

-.141250 

-.192852 

-1.4 

.162881 

.745436 

.276436 

-.097540 

-.238752 

-1.2 

.219437 

.688524 

.291034 

-.048050 

-.250453 

-1.0 

.287600 

.629686 

.295718 

.000547 

-.230966 

-  .8 

.367562 

.570849 

.291238 

.042875 

-.189622 

-  .6 

.459147 

.513695 

.279206 

.075707 

-.137889 

-  .4 

.561883 

.459534 

.261660 

.098016 

-.085824 

-  .2 

.675073 

.409261 

.240659 

.110477 

-.040233 

.0 

.797885 

.363380 

.218015 

.114769 

-.004433 

.2 

.929416 

.322069 

.195158 

.1129 46 

.020990 

.4 

1.068757 

.285262 

.173110 

.106994 

.037137 

.6 

1.215026 

.252727 

.152523 

.098587 

.045869 

.8 

1.367403 

.224132 

.133752 

.089010 

.049170 

1.0 

1.525136 

.199096 

.116935 

.079165 

.048804 

1.2 

1.687553 

.177229 

.102061 

.069643 

.046173 

1.4 

1.854058 

.158150 

.089030 

.060786 

.042308 

1.6 

2.024130 

.141506 

.077687 

.052765 

.037920 

1.8 

2.197314 

.126976 

.067859 

.045637 

.033466 

2.0 

2.373217 

.114276 

.059367 

.039387 

.029222 

2.2 

2.551498 

.103155 

.052040 

.033956 

.025338 

2.4 

2.731863 

.093396 

.045722 

.029266 

.021881 

2.6 

2.914059 

.084813 

.040273 

.025234 

.018866 

2.8 

3.097868 

.077244 

.035568 

.021774 

.016280 

3.0 

3.283101 

.070551 

.031501 

.018808 

.014096 
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